REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis 
Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a 
collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YYYY) 
31-10-2005 

2.  REPORT  TYPE 

Journal  Article 

3.  DATES  COVERED  (From  -  To) 

4.  TITLE  AND  SUBTITLE 

5a.  CONTRACT  NUMBER 

Theory  of  Monte-Carlo  Simulations  of  the  Magnetic  Circular  Dichroism  (MCD)  Spectra 
of  Alkali  Metal/Rare  Gas  Systems  (PREPRINT) 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

John  W.  Kenney  (Concordia  University);  Jerry  A.  Boatz  (AFRL/PRSP);  Heidi  A.  Terrill 

5d.  PROJECT  NUMBER 

23030423 

Vosbein  (Stennis  Space  Center) 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory  (AFMC) 

AFRL/PRSP 

10  E.  Saturn  Blvd. 

Edwards  AFB  CA  93524-7680 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AFRL-PR-ED-JA-2005-439 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S 
ACRONYM(S) 

Air  Force  Research  Laboratory  (AFMC) 

AFRL/PRS 

5  Pollux  Drive 

Edwards  AFB  CA  93524-7048 

11.  SPONSOR/MONITOR’S 
NUMBER(S) 

AFRL-PR-ED-JA-2005-439 

12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited  (AFRL-ERS-PAS-2005-277) 

13.  SUPPLEMENTARY  NOTES 

Submitted  for  publication  in  the  International  Journal  of  Quantum  Chemistry 

14.  ABSTRACT 

The  history  of  magnetic  circular  dichroism  (MCD)  spectroscopy  in  the  study  of  alkali  metal/rare  gas  (M/Rg)  cryogenic 
systems  is  reviewed  in  the  context  of  developing  a  better  understanding  of  alkali  metal/hydrogen  systems  of  current  interest  to  the 
U.S.  Air  Force  as  enhanced-performance  cryogenic  rocket  propellants.  A  new  theory  for  simulating  the  MCD  spectra  of  M/Rg 
systems  is  presented  together  with  a  careful  discussion  of  the  theory’s  implicit  and  explicit  approximations  and  their  implications. 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF  RESPONSIBLE 
PERSON 

Dr.  Scott  A.  Shackelford 

a.  REPORT 

Unclassified 

b.  ABSTRACT 

Unclassified 

c.  THIS  PAGE 

Unclassified 

A 

24 

19b.  TELEPHONE  NUMBER 

(include  area  code) 

N/A 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  239.18 


Unclassified 


Unclassified 


Unclassified 


PREPRINT 


HEP/123-qed 


Theory  of  Monte- Carlo  simulations  of  the  Magnetic  Circular 
Dichroism  (MCD)  spectra  of  alkali  metal/rare  gas  systems 

John  W.  Kenney,  III 

Chemical  Physics  Laboratory,  Concordia  University, 

1530  Concordia  West,  Irvine,  CA  92612 

Jerry  A.  Boatz 

Air  Force  Research  Laboratory,  Space  and  Missle  Propulsion  Division;  AFRL/PRSP, 
10  E.  Saturn  Blvd.,  Edwards  AFB,  CA  93524 

Heidi  A.  Terrill  Vosbein 
Naval  Research  Laboratory,  Acoustic  Simulation, 

Measurements  and  Tactics  Branch,  Code  7180, 

Bldg.  1005,  Stennis  Space  Center,  MS  39529 
(Dated:  ;  Received  text  ;  Revised  text 

Accepted  text  ;  Published  text 


This  page  is  Distribution  A:  Approved  for  public  release,  distribution  unlimited 


101 


Abstract 

The  history  of  magnetic  circular  dichroism  (MCD)  spectroscopy  in  the  study  of  alkali  metal/rare 
gas  (M/Rg)  cryogenic  systems  is  reviewed  in  the  context  of  developing  a  better  understanding  of 
alkali  metal/hydrogen  systems  of  current  interest  to  the  U.S.  Air  Force  as  enhanced-performance 
cryogenic  rocket  propellants.  A  new  theory  for  simulating  the  MCD  spectra  of  M/Rg  systems  is 
presented  together  with  a  careful  discussion  of  the  theory’s  implicit  and  explicit  approximations 
and  their  implications.  This  theory  uses  a  classical  Monte  Carlo  (MC)  simulation  scheme  to  model 
the  perturbing  effects  of  the  Rg  environment  on  the  2S  — r  2P  MCD- active  transition  of  the  M 
atom.  The  theory  sets  up  the  MC-MCD  simulation  as  a  6  x  6  matrix  eigenvalue/eigenvector 
problem  in  the  2P  manifold  in  which  is  included  the  effects  of  M-Rg  interactions,  metal  atom  spin- 
orbit  coupling  in  the  2P  manifold,  magnetic  Zeeman  perturbations  of  the  2S  and  2P  manifolds, 
Boltzmann  temperature  factors,  and  electric  dipole  transition  moment  integrals  for  left  circularly 
polarized  (LCP)  and  right  circularly  polarized  (RCP)  light.  The  theory  may  be  applied  to  any 
type  of  trapping  site  of  the  host  M  in  the  guest  Rg  matrix;  a  single  atom  substitutional  metal  atom 
trapping  site  (one  host  Rg  atom  is  replaced  by  one  guest  M  atom)  is  modeled  in  this  study  for  M 
=  Na  and  Rg  =  Ar.  Two  temperature  factors  are  used  in  these  simulations;  a  lattice  temperature 
to  model  the  mobility  of  the  Rg  lattice  and  a  magnetic  temperature  to  model  Boltzmann  factors 
in  the  2S  ground  manifold.  The  6x6  eigenvalue/eigenvector  problem  is  solved  for  a  number 
of  randomly  generated  and  suitably  averaged  Rg  configurations  to  yield  the  simulated  MC-MCD 
spectrum  for  the  single  substitutional  Na/Ar  system.  The  MC-MCD  simulations  of  Na/Ar  give  the 
characteristic  triplet  MCD  spectrum  with  the  correct  Boltzmann  temperature  dependence.  The 
simulated  MC-MCD  spectrum  correctly  inverts  when  the  direction  of  the  applied  magnetic  field 
is  reversed.  Addition  of  the  LCP  and  RCP  absorbances  gives  rise  to  a  characteristic  2S  — >  2P 
triplet  absorption  feature. 
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I.  INTRODUCTION 

Alkali  metal/rare  gas  solids  (M/Rg  solids)  constitute  an  important  class  of  model  sys¬ 
tems  of  interest  to  the  U.  S.  Air  Force  in  simulating  the  properties  of  enhanced-performance 
cryogenic  rocket  propellants  [1-4],  In  particular,  the  alkali  metal/rare  gas  (M/Rg)  sys¬ 
tems  share  many  physical,  chemical,  and  spectroscopic  characteristics  with  M/hydrogen 
systems,  where  in  a  potential  rocket  propulsion  application,  M  is  a  high  energy  density 
material  (HEDM);  e.g.,  a  low  mass  alkali  metal  or  Group  13  or  14  element  [5,  6].  Boat-z 
and  Fajardo  recently  modeled  2 S  — >  2P  electronic  absorption  spectra  in  M/Rg  systems  by 
using  quantum  mechanical  first-order  degenerate  perturbation  theory  in  conjunction  with  a 
classical  Monte  Carlo  (MC)  method  to  account  statistically  for  the  effects  of  the  Rg  per¬ 
turbations  on  the  2S  and  2P  terms  of  the  M/Rg  systems  [7].  Their  simulations,  while  not- 
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completely  quantitative,  do  recover  many  of  the  key  attributes  of  the  electronic  absorption 
spectra  of  these  systems  including  the  characteristic  three-peaked  or  "triplet"  absorption 
band  shape.  Lawrence  and  Apkarian  take  a  very  similar  approach  in  modeling  the  2P  terms 
of  halogen  atoms  doped  in  cryrogenic  matrices  (X/Rg)  [8].  Boatz  and  Fajardo  make  the 
reasonable  assumption  that  spin-orbit  coupling  represents  a  small  perturbation  of  the  2P 
term  compared  to  the  Rg  perturbation,  at  least  for  the  lighter  alkali  metals  of  interest  to 
them  as  HEDM  rocket  propellant  dopants,  and  on  this  basis  justify  the  exclusion  of  2P 
spin-orbit  coupling  terms  in  their  model.  The  Boatz  and  Fajardo  first-order  perturbation 
matrix  for  the  2P  term  is,  therefore  a  3  x  3  matrix  whose  matrix  elements  are  defined  with 
respect  to  a  suitable  zeroth-order  orbital  p  basis  in  the  angular  momentum  representation 
(| n  l  mi)  =  | n  1  1) ,| n  1  0},|n  1  —1)).  Spin  is  not  considered.  Lawrence  and  Apkarian  do 
include  a  2P  spin-orbit  term  for  the  halogen  in  their  treatment.  They  argue  convincingly 
that  the  large  halogen  2P  spin-orbit  coupling  interaction  is  not  greatly  affected  by  the  matrix 
environment,  i.e.,  the  use  of  the  atomic  spin-orbit  coupling  constant  of  the  halogen  in  the 
formalism  represents  a  good  approximation.  The  Lawrence  and  Apkarian  2P  perturbation 
matrix  is  a  6  x  6  matrix  whose  individual  matrix  elements  are  defined  with  respect  to  a 
suitable  zeroth-order  Cartesian  p  basis  with  spin  {| n  l  e  s  ms)  =  | n  l  x  \  \n  l  y  \  ^), 
I nlz\  §),  | nl  x\  -|>,  | nly\  -\),  \ nl  z  \  -|>}. 

Magnetic  circular  dichroism  (MCD)  spectroscopy,  which  measures  the  differential  ab¬ 
sorption  of  left  circularly  polarized  (LCP)  light  vs.  right  circularly  polarized  (RCP)  light 
in  a  sample  placed  in  a  magnetic  field  aligned  parallel  to  the  propagation  direction  of  the 
light,  has  a  long  and  venerable  history  of  providing  insights  to  the  nature  of  the  chemical 
and  physical  environment  surrounding  the  MCD-active  chromophore  [9].  In  particular,  the 
MCD-active  2 S  — >  2P  electronic  transitions  of  alkali  halide  F  centers  [10-14]  and  M/Rg 
systems  [15-18]  have  attracted  considerable  experimental  and  theoretical  attention  over  the 
years.  The  usual  approach  is  to  extract  parameters  such  as  the  g  value,  the  spin-orbit 
coupling  constant,  and  linear  vibrational  coupling  constants  from  a  moment  analysis  of  the 
experimental  MCD  spectrum  (or  spectra)  [9,  11].  Only  a  few  full  MCD  spectral  simulations 
exist  in  the  literature.  These  simulations  of  2S  —*■  2P  MCD  spectra  rely  upon  dynamic 
Jahn- Teller  models  in  which  an  a  priori  choice  is  made  as  to  which  specific  interaction  vi¬ 
brational  modes  will  act  to  lift  the  degeneracy  of  the  2P  term  [19,  20].  This  simulation 
approach  is  of  limited  value  in  modeling  the  MCD  spectra  of  M  atoms  trapped  in  novel 
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and  heretofore  unobserved  sites  in  Rg  matrices  for  which  specific  site  geometries  and  lattice 
mode  types  are  as  yet  unknown.  A  MC  approach  for  handling  Rg  perturbations  in  MCD 
simulations  of  M/Rg  systems  naturally  suggests  itself.  A  full  development  of  the  theory 
of  Monte  Carlo-MCD  (MC-MCD)  simulations  of  2  S'  — >  2P  MCD  spectra  is  presented  in 
the  sections  that  follow.  The  MC-MCD  simulation  method  is  a  first  order  degenerate  per¬ 
turbation  method  that  straightforwardly  extends  the  MC  absorption-  simulation  methods 
of  Boatz  and  Fajardo  [7]  and  Lawrence  and  Apkarian  [8]  by  (1)  including  both  spin-orbit 
and  Zeeman  perturbations  in  the  Hamiltonian  and  (2)  separately  treating  the  2S  —>■  2P 
transitions  associated  with  LCP  and  RCP  light.  The  MC-MCD  perturbation  matrix  for 
the  2P  term  is  a  6  x  6  matrix  whose  individual  matrix  elements  are  expressed  in  terms  of 
a  suitable  uncoupled  angular  momentum  p  basis  set  with  spin  {| nlmisms)  =  |nl  1  \ 

|wl0H>»  \nl  -1  H)>  \nll  \  -|>»  |wl0§  -§}>  \nl  -!  I  -Di¬ 

ll.  THEORY 

As  noted  above,  the  theory  of  MC-MCD  spectroscopy  of  2S  — >  2P  excitations  in  M/Rg 
systems  represents  an  extension  of  the  formalism  for  MC-absorption  spectroscopy  developed 
by  Boatz  and  Fajardo  [7]  (for  M/Rg  systems)  and  Lawrence  and  Apkarian  [8]  (for  X/Rg 
systems).  Key  aspects  of  both  the  Boatz  and  Fajardo  and  Lawrence  and  Apkarian  notat-ional 
schemes  are  adopted  in  the  following  formal  development. 

A.  Total  Hamiltonian  and  zeroth-order  eigenvectors/eigenvalues 

The  total  Hamiltonian  for  an  alkali  metal  M  atom  imbedded  in  a  cryogenic  rare  gas 
matrix  comprised  of  N  Rg  atoms  in  the  presence  of  an  applied  external  magnetic  field  is 
given  by 

H  =  HM(r )  +  Hint( r;  R1;  R2, ...,  Rjv)j  (1) 

where 

Hm( r)  =  Tm{y)  +  Vm(?)  (2) 

is  the  one-electron  Hamiltonian  for  the  optically  active  electron  on  M  at  position  r,  comprised 
of  a  kinetic  energy  component  Tm( r)  and  an  electrostatic  potential  energy  component  Vjw(r). 
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The  potential  14/  (r )  is  the  potential  experienced  by  the  optically  active  electron  in  the  free 
gas-phase  M  atom  in  the  absence  of  spin-orbit  coupling  or  an  external  magnetic  field. 

To  describe  completely  the  2  S'  ground  term  of  M,  two  zeroth-order  eigenvectors 
{ |  a  0  0  |  ms ) ;  vns  =  1/2.  -1/2)  are  needed;  for  the  2P  excited  term,  six  zeroth-order  eigen¬ 
vectors  {| nlmisms);  mi.  =  1,  0,  -1,  ms  =  1/2,  -1/2}  are  required.  An  uncoupled  angular 
momentum  representation  {\nl  mi  sms)}  for  the  zeroth-order  eigenvectors  is  chosen,  with¬ 
out  loss  of  generality,  because  of  the  transparent  way  in  which  MCD  selection  rules  may  be 
expressed  in  terms  of  this  representation  ( vide  infra).  However,  it  is  equally  valid  to  use 
zeroth-order  eigenvector  sets  expressed  in  the  coupled  angular  momentum  representation 
{| nlsjm)  or  the  Cartesian  representation  {| nlesms)}  (e  =  x,y,z  for  l  =  1)  since  the 
uncoupled,  coupled,  and  Cartesian  representations  are  mutually  related  by  unitary  trans¬ 
formations  [21].  In  this  treatment,  the  interaction  Hamiltonian  is  comprised  of  four  terms, 
each  of  which  depends  upon  the  nuclear  position  vectors  {74}of  all  of  the  N  Rg  atoms  in  the 
matrix;  a  metal-rare  gas  interaction  term,  a  rare  gas-rare  gas  interaction  term,  a  spin-orbit 
term,  and  a  Zeeman  term, 

Hint (r;  Ri,  R2, ...,  R/v)  =  14/-i?9(r;  Ri,  R2,  •••,  Rat)  +  VRg-Rg(r'i  Ri,  R2,  ••••  R n)  (3) 

T  Hso(Tj  Ri,  R2)  ■  ■  ■  ■  R/v)  T  HZeeman{,r,  R"  1  ■  R2j  ■■■)  Rjv) 

The  M  nucleus  is  taken  to  be  the  origin  of  the  coordinate  system  with  respect  to  which  r  and 
the  {-Rfc}  are  defined.  In  this  treatment,  it  is  assumed,  following  Lawrence  and  Apkarian, 
that  M-Rg  interactions  do  not  affect  Rg-Rg  interactions  [8];  no  modifications  have  been 
made  to  the  Rg-Rg  potentials  to  account  for  perturbations  arising  from  the  presence  of  the 
M  atom. 

B.  The  Rg-Rg  and  M/Rg  potentials 

The  Rg-Rg  interaction  potential  in  Eq.  3  can  be  approximated  as  a  simple  sum  over  the 
appropriate  S  symmetry  Rg-Rg  diatomic  potentials 

VRg-Rg{ Rl,  R2,  •••,  R/v)  =  V Rg-Rg ( | Rm  —  Rfc|)  =  Ur9-R9  (4) 

m=l,N  m<k 

following  the  notation  of  Boatz  and  Fajardo  [7].  It  should  be  noted  that  URg-Rg  has  the 
same  value  for  both  the  2 S  and  2P  terms  at  this  level  of  approximation. 
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The  M-Rg  potential,  which  in  general  exhibits  a  complicated  dependence  upon  the  Rg 
positions  (Rk),  is  expressed  in  terms  of  a  Legendre  polynomial  expansion  as 

VM-Rg(r 5  Rl,  R-2,  •••,  Rjv)  =  E  E  VL(r,  Rk)PL(9, 0;  $k,  fa)  (5) 

k=l,N  L=0, infinity 

where  r  =  (r,  9,  fa)  is  the  optically  active  electron  position  vector,  Rk  =  (Rk:9k,fa)  is  the 
position  vector  for  Rg  nucleus  k.  Vfar,  R*,)  is  the  radial  potential  function  of  order  L,  and 
Pl(0,  fa  dk,  4>k)  is  the  Legendre  polynomial  of  order  L 

Pl(0,  fa  0k,  fa)  =  4tt(2 L  +  l)"1  ^  Ylm(6,  <i>)Y£M(0k,  fa)  (6) 

M=—L,L 

and  the  Ylm  arc  spherical  harmonic  functions  expressed  in  the  standard  Condon  and  Shortley 
phase  convention  [22],  Eq.  6  is  an  expression  of  one  of  the  possible  forms  of  the  spherical 
harmonic  addition  theorem  [21]. 

C.  Spin-orbit  coupling  operator 

The  spin-orbit  coupling  operator  for  the  optically  active  electron  may  be  written,  using 
the  standard  defining  equation  for  the  spin-orbit  interaction  [23],  in  terms  of  the  potential 
gradient  of  the  free  M  and  M-Rg  interaction  potentials 

Hso(y- R1;  R2,  ...,  Rjv)  =  HSO)M(y )  +  Hso,M-Rg(r'i  Rl,  R-2)  •••,  Rat) 

=  {(2me2c2r)_1  V[VM(r)  +  VM-Rg(r:  Ri,  R2, Rat)]  x  p}  ■  S 
=  Kiif(r)  +  £,M-Rg(r)]L  ■  S  +  [(2me2c2r)_1  (7) 

xv(E  E  VL(r,Rk)PL(6,<l>-,0k,fa)}xp]-S 

k=l,N  L=l, infinity 

where  p  and  S  are,  respectively,  the  linear  momentum  and  spin  operators  for  the  optically 
active  electron,  L  =  r  x  p  is  the  angular  momentum  operator  for  the  optically  active  electron 
and 

€m(t)  =  (2me2c2r)~1d[VM(r)]/dr  (8) 

is  the  standard  spin-orbit  coupling  operator  for  the  free  M  atom.  The  M-Rg  interaction 
potential  of  Eq.  5  can  be  broken  up  into  spherical  and  non-spherical  components 

Ym~ Rg (r ;  Ri,  R2,  •••,  Rat)  =  vo(r’Rfc)+  S  FL(r,R k)PL(0,<f>]0k,fa).  (9) 

k=l,N  k=l,N  L=l, infinity 
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In  the  form  of  the  theory  currently  in  use  for  MC-MCD  simulations,  it  is  assumed  that 
t-he  gradient  expression  in  Eq.  7  is  small  and  can  be  neglected;  i.e.,  no  L  >  0  terms  are 
assumed  to  contribute  to  the  spin-orbit  coupling  interaction.  The  total  spin-orbit  coupling 
operator  of  the  M/Rg  system  in  the  spherical  approximation  is  thus  given  by 

Hso{ r)  =  HsoM( r)  +  H30,m-r9( r)  =  [U(r)  +  Cm-rM L  •  S  =£(r,  L  •  S)  (10) 

where  the  dependence  on  the  Rg  position  vectors  {Rk},  while  present  is  not  explicitly  shown 
in  the  interest  of  preserving  a  compact  notation  in  the  final  form  of  the  equations  used. 

No  spin-orbit  coupling  is  possible  when  the  optically  active  electron  is  in  the  2  S'  ground 
state  manifold  where  1  =  0.  In  the  2P  excited  state  manifold  where  1  =  1,  the  spin-orbit 
coupling  constant  is  expressed  formally  as  the  sum  of  the  radial  expectation  values  of  the 
spherically  symmetric  M  and  M-Rg  spin-orbit  coupling  operators 

Cm  =  (n  1 1  £(r)  | n  l)  =  (n  l\ £M{r)  \n  l)  +  (n  l\ CM-Rg(r )  \ n  l)  =  Cm  +  Cm-rg  (H) 

If  the  radial  eigenvector  | n  l)  is  the  free  atom  P  state  atomic  radial  eigenvector  for  M,  or  if 
| n  l)  closely  approximates  this  eigenvector,  then  Cm  is  the  free  alkali  atom  spin-orbit  coupling 
constant  and  Cm-rg  is  the  spherical  perturbation  to  the  spin-orbit  coupling  constant  induced 
in  M  by  the  Rg  matrix  environment.  Eq.  11  provides  an  insightful  context  for  discussing  the 
possibilities  of  (1)  spin-orbit  coupling  constant  sign  reversal  or  (2)  reduction  of  the  magnitude 
of  a  spin-orbit  coupling  constant  (partial  quenching,  in  a  2P  M-Rg  system.  For  any  given 
alkali  metal  M  the  free  atom  2P  spin-orbit  coupling  constant  is  positive;  i.e.,  Cm  >  0.  Thus, 
for  a  2P  spin-orbit  coupling  constant  Cm  to  be  negative  in  a  M-Rg  system,  the  M-Rg  spin 
orbit  contribution  must  be  negative  and  larger  in  magnitude  than  the  free  M  atom  spin- 
orbit  coupling  constant,  i.e.,  Cm-rg  <  0  and  \Cm-rg\  >  Cm  The  conditions  for  the  partial 
quenching  case,  where  Cm  is  still  positive  but  smaller  than  Cm  are:  Cm-rg  <  0  and  |  Cm-rg  \  < 
Cm-  In  the  actual  MC-MCD  simulations,  Cm  is  introduced  as  an  adjustable  parameter.  The 
two  cases  may  be  explored  in  the  MC-MCD  simulations  by  making  appropriate  choices  for 
Cm  with  respect  to  the  known  positive  Cm  of  the  alkali  metal  M  atom  of  interest. 

D.  Zeeman  Perturbation 

The  Zeeman  perturbation  operator  is  given  by 

H Zeemaniy,  R-l,  R2,  Rat)  =  iocal  (L  +  2S)  (12) 
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where  // «  is  the  Bohr  magneton,  L  and  S  are  the  orbital  and  spin  angular  momentum 
operators  for  the  optically  active  electron,  and  B iocai  is  the  magnetic  field  locally  experienced 
by  the  optically  active  electron  as  a  result  of  the  application  of  an  external  magnetic  field 
B  to  the  M/Rg  system.  Both  the  magnitude  and  the  direction  of  the  local  magnetic  field 
vector  are  influenced  by  the  positions  of  the  Rg  atoms  and  by  the  nature  of  the  electronic 
state  of  the  optically  active  electron  (i.e. ,  2  S'  or  2P).  The  relationship  between  the  external 
magnetic  field  and  the  local  magnetic  field  thus  may  be  expressed  as 

B;ocoz  =  B/ocaz(Ri,  R2, ....  Rjv)  =  [1  —  crn/(Ri,  R2, ...,  Rjv)]B  (13) 

where  <rnj(R1,R2,  ...,Rjv)  is  a  3  x  3  shielding  tensor  whose  individual  components  depend 
upon  both  the  nature  of  the  electronic  state  being  shielded  and  the  positions  of  the  Rg  atoms 
in  the  matrix.  In  actual  MC-MCD  simulations,  it  is  assumed  that  the  local  magnetic  field 
is  aligned  along  the  laboratory  z  axis 


®  local  Biocal^z  (14) 

which  gives 

Hzeeman  /i B B  (1  (Tni )  (T2  T  2 S*2)  (15) 

Different  values  of  Biocai  can  be  chosen,  if  desired,  to  model  the  effects  of  slightly  different 
2S  and  2P  shielding  constants  an0  and  anl.  Changes  in  the  direction  of  Biocai  away  from  the 
z  axis  arising  from  the  tensor  character  of  ani  can  be  modeled  by  preserving  the  laboratory 
z  axis  orientation  of  the  Zeeman  perturbation,  which  is  very  desirable  from  the  standpoint 
of  matrix  element  evaluation,  and  rotating  the  Rg  position  vector  set  through  an  arbitrary 
unitary  rotation  transformation  U { R/, }  =  {R/}.  In  any  event,  it  is  expected  that  differences 
between  B  and  B/oca/  will  be  comparatively  small  for  closed  shell  Rg  atoms  in  their  ground  1S 
electronic  states.  Certainly  within  the  chosen  MC-MCD  simulation  restrictions  mentioned 
earlier  of  neglecting  polarization  effects  and  the  higher  excited  states  of  the  Rg  atoms,  B 
will  not  differ  much  from  Biocai. 

E.  Diagonalization  of  the  full  Hamiltonian  matrix  in  the  2S  and  2P  manifolds 

An  implementation  of  degenerate  first-order  perturbation  theory  is  required  to  arrive  at 
approximate  2 S  and  2P  eigenvalues  and  eigenvectors  of  the  full  M/Rg  Hamiltonian  given  in 
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Eq.  1.  These  eigenvalues  and  eigenvectors  contain  the  effects  of  the  various  terms  in  the 
interaction  Hamiltonian  Eq.  3.  First-order  perturbation  theory  gives  eigenvalues  that  are 
correct  through  first  order.  The  associated  eigenvectors  are  those  particular  zeroth-order 
eigenvectors,  expressed  as  orthonormal  linear  combinations  of  the  zeroth-order  eigenvector 
set  {\nlmi  sms)}, 

dh  =  \nli)  =  Cfm,am'\nlmisme),  (16) 

mi=—l,l  ms=— l/2,l/2 

(nli'\nli)  =  5(i' ,i),  (17) 

that  diagonalize  the  2x2  ( 2  S'  manifold)  or  6  x  6  (2P  manifold)  matrix  formed  by  sand¬ 
wiching  the  full  Hamiltonian,  Eq.  1,  within  the  same  uncoupled  zeroth-order  eigenvector  set 
{| nlmisms)};  i.e., 


Both  the  first-order  eigenvalues  and  the  associated  zeroth-order  eigenvectors  are  needed 
for  MC-MCD  simulations.  Since  the  first  matrix  element  on  the  right  side  of  Eq.  18  is  the 
zeroth-order  matrix  element  of  the  free  M  Hamiltonian,  it  is  automatically  diagonal  in  either 
the  original  zeroth-order  basis  {  \nlmisms)  } 


(n'l'm!l  s'm!s\HM\n Imi  sms)  =  £niS(n' ,n)<5(//,/)<5(m(,m;)5(s/,  s)5(m's,ms)  (19) 

or  the  new  zeroth-order  basis  { |n  li)  } 


( nli'\HM\nli )  =  £ni5(n'  ,n)5(V  ,l)6(ir  ,i).  (20) 

Within  the  2S  manifold  the  interaction  matrix  on  the  right  side  of  Eq.  18  also  is  diagonal, 
thus 


T,,(2P)  =  |nO*>=  ^  C\ 


n00ms 


ms=- 1/2, 1/2 


n  0  0  -  ms 
2  , 


(21) 


which  reduces  to 


^i(2S) 

^2(2^) 


1  1\ 

n  °  °  2  2  /  ’ 

1  1  \ 
n  0  0  -  —  )  . 
2  2/ 


(22a) 

(22b) 
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The  i  index  is  mapped  onto  the  ms  index  one-to-one  in  this  specific  case.  The  final  expres¬ 
sions  for  the  2S  eigenvalues,  correct  through  first  order,  are  obtained  by  collecting  together 
the  matrix  elements  previously  derived  for  the  Rg-Rg,  M-Rg,  and  Zeeman  perturbations: 


Ei(2S)=sn0CS)  +  En0iCS ) 

=  ZnoCS)  +  U Rg-RgC  S)  +  U M-RgC  S)  +  UzeemanC  $) 

=  £noCS)  +  E  E  y Rg-Rg, XT,  (|Rm-Rfc|)  (23) 

m=l,N  m<k 

+  E  E  VM-Rg,XY,{Rk)  +  HbB  localisms). 

m=l,N  m<k 

In  Eq.  23  the  eigenvalue  convention  of  Boatz  and  Fajardo  [7]  is  adopted:  Er  is  an  eigenvalue 
of  a  specific  2 S  state  T,  correct  through  first  order,  en0  is  the  zeroth  order  2S  eigenvalue 
for  the  optically  active  electron  of  the  free  M  atom,  and  En0i  is  the  first-order  eigenvalue 
correction  term  arising  from  the  interaction  Hamiltonian. 

To  obtain  2P  final  state  eigenvalues  correct  through  first  order 

Ef(2P)  =  £„1  (2P)  +  Em,(2P)  (24) 


the  first-order  eigenvalue  set  ( Enif(2P );  f  =  1,6)  must  be  calculated  by  diagonalizing  the  6  x 
6  2P  interaction  matrix  in  Eq.  18.  The  individual  matrix  elements  of  this  interaction  matrix, 
shown  below  in  Table  1,  are  obtained  by  collecting  together  the  matrix  elements  previously 
derived  for  the  2P  Rg-Rg,  M-Rg,  spin-orbit,  and  Zeeman  perturbations.  It  should  be  noted 
that  the  Pug-iig  and  the  (Vq) M_Rg  perturbations  are  diagonal  perturbations  in  this  6x6 
matrix  and  can  be  separated  from  the  6x6  matrix  eigenvalue  problem  as  illustrated  by 
Boatz  and  Fajardo  in  their  Eqs.  (6)  and  (8)  [7].  The  associated  final  state  2P  eigenvectors 
that  diagonalize  this  matrix,  are 


mI?P)  =  \nlf)=  Y,  E 

rai=— 1,1  ms=- 1/2, 1/2 


c 


nlm\  I 


/ 


U  1  777-1  -  TTls 


(25) 


where  /  =  1,6  and 


(nlf\nlf)  =  <5(/',  /). 


(26) 
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F.  2S  — >  2P  transition  energies  and  transition  moment  integrals  for  M/Rg  MCD 
spectroscopy 

The  combined  effects  of  the  Rg-Rg,  M-Rg,  spin-orbit,  and  Zeeman  perturbations  on  the 
optically  active  electron  of  a  M/Rg  system  act  to  split  the  ground  2S  manifold  into  two 
energetically  distinct  states  {T,;  i  =  1,2}  and  the  excited  2P  manifold  into  six  energetically 
distinct  states  { VR / ;  /  =  1,6}.  The  specific  characteristics  of  the  2S  and  2P  eigenval¬ 
ues  and  their  associated  eigenvectors  are  dependent  upon  the  particular  Rg  configuration 
{Ri,  R2, Rjv}  within  which  the  problem  is  solved.  Thus,  for  a  given  Rg  configuration 
{Ri,  R2, ...,  Rjv}  and  magnitude  of  the  magnetic  field,  up  to  12  energetically  and  optically 
distinct  transitions  may  be  identified  within  the  manifold-to-manifold  2 S  — >  2P  transition 
of  the  M/Rg  system.  The  transition  energies  of  these  12  2S  — >  2P  transitions  are  given  as 

hv„  =  E,CP)  -  Ei(2S)  (27) 

where  i  =  1,2  and  /  =  1,6.  The  Rg-Rg  interaction  energy  terms  cancel  in  Eq.  27  since 

U Rg-Rg  (2P)  =  IJ He/  Rq  ( 2  S )  =  URg_Rg. 

Electronic  transitions  involving  the  optically  active  electron  of  a  M/Rg  system  are  well 
described  as  electric  dipole  transitions.  In  MCD  spectroscopy,  each  of  the  12  transitions  may 
be  induced  to  a  greater  or  lesser  extent  by  either  LCP  or  RCP  light  interacting  with  the 
M/Rg  system  in  the  electric  dipole  limit-  and  tuned  to  the  particular  transition  energy  }wf%. 
Thus,  both  LCP  and  RCP  electric  dipole  transition  moment  integrals  must  be  computed  for 
each  of  the  12  possible  combinations  of  initial  and  final  states,  i.e.,  24  transition  moment 
integrals,  12  for  LCP  light  and  12  for  RCP  light,  are  needed  to  describe  the  MCD  spectrum 
of  a  25'  — >  2P  transition  in  M/Rg  system  in  a  particular  Rg  configuration  {Ri,  R2, ...,  Ry}; 
i.e., 

<'P/(2P)|OLCp|<I',(2S))  (28a) 

(®/(2P)|Obcp|Ki(2S))  (28b) 

where  i  =  1,2  and  /  =  1,6.  Olcp  and  0Rcp  are  the  electric  dipole  transition  operators 
for  LCP  and  RCP  light,  in  Eqs.  28a  and  28b.  The  electric  dipole  moment  operator  m 
of  the  M/Rg  system  is,  to  a  very  good  initial  approximation,  determined  by  the  position 
r  =  ( x ,  y,  z )  and  charge  — e  of  the  optically  active  electron.  The  Rg  atoms  are  electrically 
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neutral  and  their  contributions  to  the  dipole  moment  can  be  neglected  if  polarization  effects 
are  assumed  to  be  small.  The  M+  core  is  positioned  at  the  origin  of  the  chosen  coordinate 
system  and  thus  does  not  contribute  to  the  expression  for  the  dipole  moment  at  this  level 
of  approximation,  i.e., 

m  =  ^  qin  =  -er  +  eO  =  -e(xex  +  yey  +  zez )  +  e(0ex  +  Oe,,  +  Oe.) 

i 

=  —e(xex  +  yey  +  zez).  (29) 

However,  the  optically  active  electron  at  position  r  will  induce  a  slight  polarization  in 
the  M+  core  and  the  Rg  matrix.  This  polarization  can  be  represented  as  a  small  negative 
electrostatic  image  charge  —q'  at  — r.  Thus,  polarization  effectively  reduces  the  charge  of 
the  optically  active  electron  to 

—e'  =  —  e  +  q'  (30) 

where  |  —  e'\  <  |  —  e|.  and  the  dipole  moment  operator  of  the  M/Rg  system  becomes 

m  =  — e'r  =  —e'(xex  +  yey  +  zez).  (31) 

While  it  can  be  argued  that  the  image  charge  —q'  (and  hence  —  e')  is  not,  strictly  speaking, 
a  constant,  but  will  exhibit  both  a  radial  dependence  (i.e.,  different  q'  for  small  r  vs.  large  r) 
and  angular  dependence  arising  from  the  specific  positions  of  the  Rg  atoms  (Ri,  R2, ...,  Rjv) 
in  the  M/Rg  system,  these  effects  are  neglected  at  this  level  of  the  theory.  Within  the  Piepho 


and  Schatz  unit  vector  conventions  for  LCP  and  RCP  light  [9] 

ttlcp  =  2_1/2(ea,  -  iey )  (32a) 

nRcp  =  2_1/2(ex  +  iey )  (32b) 

the  LCP  and  RCP  electric  dipole  moment  operators  for  the  M/Rg  system  become 

OLCp  =  m  •  n*LCP  =  -e'2 ~1/2(x  +  iy)  =  e' (An ,4>)r,  (33a) 

0RCp  =  m  •  n*RCP  =  -e'2-1/2(x  -  iy )  =  -e' (An /3)1/2yi_i(6>,  <j>)r.  (33b) 


Since  each  of  the  T,(2,S')  and  f(2P)  states  is  expressed  in  terms  of  the  uncoupled  zeroth- 
order  angular  momentum  eigenvector  set  {| nlmisms)}  (see  Eqs.  16  and  25),  the  specific  form 
of  the  transition  moment  integrals  given  in  Eqs.  33a  and  33b  may  be  the  determined  by 
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evaluating  integrals 


(n  l'm[  s'w!s\OLCp\nlmi  sms)  =  ( nl'm'll/2m,s\OLCp\n001/2ms )  =  /(5(m),  l)^???/,  m,,), 

(34a) 

(nl'm!l  s'm's\0Rcp\nlmisms)  =  ( nl,m,ll/2m,s\Opcp\n^01/2ms )  =  —  f5(m'u  —  l)5(m's,  ms) 

(34b) 

where 

Jlcp  =  f  =  (e' [4tt / 3] 1/2)  (nl|r|nO)  (Yn|Yn|Y0o)  =  (e'[3]”1/2)  (nl|r|nO)  (35a) 

Ircp  =  ~f  =  (— e' [4tt/3]1/2)  (nl|r|nO)  (Yi_i|Yi_i|Y0o)  =  (e'[3]_1/2)  (nl|r|nO)  (35b) 

Eqs.  34a  and  35a  express  the  fact  that  an  LCP  electric  dipole  transition  acts  on  the  orbital 
angular  momentum  component  of  the  uncoupled  wave  function,  increasing  both  l  and  mi 
by  one  unit  while  leaving  s  and  ms  unchanged.  The  corresponding  RCP  electric  dipole 
transition  described  in  Eqs.  34b  and  35b  also  acts  on  the  orbital  angular  momentum  in  the 
uncoupled  wave  function,  increasing  l  by  one  unit  and  decreasing  to/  by  one  unit.  In  Eqs. 
35a  and  35b  /  is  the  magnitude  of  the  evaluated  transition  moment  integral.  Note  t-hat- 
the  value  of  /  is  the  same  for  all  allowed  transitions.  By  combining  the  transition  moment 
integral  expressions  in  Eqs.  28,  34,  and  35  with  the  appropriate  expressions  for  T,(2S')  and 
X\>f(-P)  in  the  uncoupled  basis  as  given  in  Eqs.  16,.  17,  and  25,  final  expressions  for  the 
contribution  to  the  MCD  spectrum  by  the  transition  i  — >  /  at  energy  fiuj,  can  be  written: 

AA  '(vfi)  =  ALCp{pfi )  —  ARCp(vfi) 

= a{  {yfep)\oLCPm2s))  r  - 1  <>m2p)|oRCp|<»,(2s))  n 

=  a/2{|C},1'm‘"1’m-|2  -  1C}*1™— 1-"*-|2)  x  [1  +  (e-2<‘»iJ/‘T  -  1  )S(m„  1/2)]  (36) 

where  ck  is  a  numerical  constant  that  relates  the  squared  transition  moment  integrals  to  the 
absorbances  and  the  initial  state  index  i  is  associated  with  ms  =  —1/2  or  1/2  in  the  2S  mani¬ 
fold.  The  MCD  spectral  amplitude  at  energy  h,Ufl  for  the  Rg  configuration  {Ri,  R2, ...,  R/v} 
is  seen  to  be  related  to  differences  of  squares  of  selected  2P  eigenvector  coefficients  com¬ 
puted  for  that  configuration.  The  term  in  square  brackets  on  the  right  side  of  Eq.  36  is  a 
Boltzmann  factor  for  the  2S  manifold.  The  Zeeman  perturbation  splits  the  energies  of  the 
ms  =  —1/2  and  1/2  states  of  the  2 S'  manifold  by  a  factor  of  2///Ri.  The  ms  =  —1/2  and 
1/2  states  will,  therefore,  have  different  thermal  populations.  This  must  be  accounted  for 
in  computing  the  various  contributions  to  the  MCD  spectrum. 
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III.  MONTE  CARLO  SIMULATION  OF  THE  MCD  SPECTRUM 


The  full  MCD  spectrum  for  a  2S'  — >  2P  transition  in  a  M/Rg  system  in  a  specific  Rg 
configuration  {Ri,R2,  R;y}  is,  therefore,  a  "stick  MCD  spectrum"  consisting  of  12  lines, 
appropriately  placed  on  the  x  (energy)  axis  at  transition  energies  hisf,.  whose  magnitudes 
and  directions  (+  or  -)  on  the  y  (AA/)  axis  are  determined  by  computing  Eq.  36  for  all 
possible  choices  of  i  and  /.  To  simulate  a  real  MCD  spectrum,  many  stick  MCD  spectra, 
each  arising  from  a  different  Rg  configuration,  are  averaged  together  using  a  simple  binning 
(histogram)  technique  based  upon  the  original  Metropolis  et  al.  algorithm  [24]  exactly  as 
implemented  by  Boatz  and  Fajardo  for  the  case  of  electronic  absorption  spectroscopy  [7].  A 
MC-MCD  simulation  begins  by  choosing  an  initial  Rg  configuration  {R1;  R2, ..., initial, 
typically  ,  representing  an  idealized  trapping  site  (e.g.,  a  single  substitutional  site  with  O/, 
symmetry).  The  MC  energy  optimization  scheme  operates  in  the  2 S'  ground  state  manifold. 
At  this  stage,  the  very  small  Zeeman  perturbation  of  the  2S  manifold  is  neglected.  MCD  stick 
spectra  arising  from  more  favorable  configurations  are  weighted  more  heavily  in  the  averaging 
process  than  those  arising  from  the  less  favorable  configurations.  In  actual  practice,  the  LCP 
and  RCP  contributions  to  the  MCD  spectrum  appearing  in  Eq.  36  are  accumulated  and 
stored  separately.  This  allows  the  electronic  absorption  spectrum  to  be  recovered  from  the 
MC-MCD  simulation  as  [9] 


^  —  A  —  ^  [A'lcp  +  A'rcpI  •  (37) 

Eq.  37  expresses  the  fact  that  the  electronic  absorption  spectra  of  M/Rg  systems  are  essen¬ 
tially  unaffected  by  the  application  of  an  external  magnetic  field. 

A.  Orientational  averaging  in  MC-MCD  simulations 

In  an  MCD  experiment,  the  propagation  direction  of  the  light,  which  is  always  parallel 
to  the  B  vector,  rigorously  defines  the  laboratory  frame  z  axis.  The  laboratory  frame  x  and 
y  axes  are  set  parallel  to  the  sapphire  deposition  window  and  perpendicular  to  B  with  the 
origin  at  the  M  nucleus.  All  Rg  positions,  eigenvectors,  and  selection  rules  are  described  with 
respect  to  this  laboratory  frame.  To  model  the  effects  of  arbitrary  Rg  lattice  orientation 
relative  to  the  laboratory  frame,  the  6x6  2P  matrix  listed  in  Table  1  is  separated  into 
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M/Rg  and  spin-orbit  +  Zeeman  matrices 

[M/Rg  +  SO  +  Zeeman]  =  [M/Rg]  +  [S'O  +  Zeeman],  (38) 

Arbitrary  orientations  of  the  Rg  matrix  relative  to  the  laboratory  frame  can  be  modeled  by 
re-writing  Eq.  38  as 

[M/Rg  +  SO  +  Zeeman]  =  [U]  [M/Rg]  [U]^1  +  [SO  +  Zeeman]  (39) 

where  U  is  an  arbitrary  randomly  generated  unitary  rotation  matrix  (e.g.,  an  Eulerian 
rotation  matrix).  The  MC-MCD  simulation  process  is  carried  out  as  before  using  Eq.  39 
rather  than  Eq.  38. 

B.  Lattice  temperature  and  magnetic  temperature  in  MC  MCD  simulations 

Temperature  effects  appear  in  two  distinct  places  in  MC-MCD  simulations.  A  lattice 
temperature  must  be  specified  for  the  MC  simulation  process.  It  may  be  argued  that  a 
proper  representation  of  quantum  mechanical  zero  point  vibrational  motions  of  the  M/Rg 
lattice  at  temperature  T  requires  the  classical  MC  simulation  to  be  carried  out  at  a  classical 
lattice  temperature  T'  >  T  [8].  However,  a  temperature  for  the  Boltzmann  factor  of  the 
Zeeman-split-  2 S  term  also  must  be  included.  This  magnetic  temperature,  Tmag,  always 
should  be  the  actual  system  temperature  T  (Tmag  =  T ).  The  MC-MCD  simulation  code 
has  the  option  of  specifying  separate  values  for  the  lattice  temperature  V  and  the  magnetic 
temperature  Tmag.  Lawrence  and  Apkarian  discuss  the  relationship  between  T'  and  T  in 
classical  MC  simulations  of  doped  Rg  systems  [8].  In  practice  T'  is  calibrated  to  T  by 
comparing  the  Rg  MC  equilibrium  lattice  spacing  (for  first  and  second  nearest  neighbors) 
with  the  experimentally  measured  values  [25].  For  Ar  at  10  K,  the  first  nearest  neighbor 
analysis  predicts  a  lattice  temperature  of  31 K,  while  the  second  nearest  neighbor  analysis 
predicts  a  slightly  lower  lattice  temperature  of  27 K.  Consequently  a  simulation  lattice 
temperature  in  the  vicinity  of  30/7  represents  the  best  modeling  of  the  quantum  zero  point- 
vibrations  of  the  Ar  lattice.  The  Na/Ar  system  at-  a  magnetic  temperature,  Tmag  =  10 A.  a 
Na  2P  spin-orbit  parameter  of  17  cm-1  (the  value  for  unperturbed  atomic  Na)  is  illustrated 
in  Figures  la  and  lb  in  which  are  presented  the  simulated  absorption  and  MCD  spectra, 
respectively,  for  three  different  choices  of  lattice  temperature,  T'  =  10/7.  30 K  and  45/7. 
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While  the  characteristic  3-peaked  "triplet"  absorption  feature  and  the  corresponding  down, 
down,  up  MCD  triplet  feature  are  preserved  for  all  three  choices  of  T'\  the  spectral  breadth 
increases  as  T'  increases,  improving  agreement  between  simulated  and  experimental  spectra. 
The  difference  in  peak  height  is  due  to  the  fact  that  the  total  area  under  the  absorption 
curve  is  normalized  to  1.0;  broader  peaks  will  be  shorter. 

C.  Other  important  properties  of  MC-MCD  simulations 

For  the  single  substitutional  site  in  the  Na/Ar  system,  absorption  and  MC-MCD  simu¬ 
lations  correctly  reproduce  other  observed  properties  of  experimental  absorption  and  MCD 
spectra  M/Rg  systems.  For  example,  the  observed  insensitivity  of  the  experimental  2S  — >  2P 
M/Rg  triplet  absorption  feature  to  the  magnetic  field  strength  and  Tmag  is  faithfully  repro¬ 
duced  in  the  MC-MCD  simulation  as  shown  in  Figure  2.  Reversal  of  the  magnetic  field 
direction  (B  =  Bez  changed  to  B  =  — Bez )  in  an  MCD  experiment  has  the  effect  of  invert¬ 
ing  the  MCD  spectrum.  If  the  normal  MCD  spectrum  follows  a  down-down-up  pattern, 
the  MCD  spectrum  with  the  reversed  magnetic  field  will  be  the  mirror  image  of  the  original 
spectrum  with  an  up-up-down  pattern.  In  Figure  3,  it  is  seen  that  the  MC-MCD  simula¬ 
tions  of  Na/Ar  correctly  gives  an  inverted  spectrum  when  B  is  replaced  with  — B.  These 
successes  suggest  the  MC-MCD  simulation  embodies  much  of  the  essential  physics  of  the 
Na/Ar  system  and  its  interaction  with  external  magnetic  fields,  the  perturbing  effects  of 
the  Rg  matrix  and  electromagnetic  radiation.  Overall,  the  MC-MCD  simulated  Na/Ar 
absorption  and  MCD  spectra  show  good  qualitative  agreement  with  experimental  Na/Ar 
spectra. 
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Fig.  1.  Comparison  of  the  temperature  effects  in  the  (a)  absorption  and  (b)  MCD  spectra 
for  the  case  of  SO  =  17 cm-1,  B  =  1.0  T .  magnetic  temperature,  Tm  =  10  K.  at  simulation 
temperatures  of  10,  30,  and  45  K 

Fig.  2.  Absorption  spectra  with  (a)  normal  and  (b)  reversed  B-fields.  1-atom  vacancy, 
SO  =  17  cm^1,  and  Tmag  =  6,  8,  10  K.  The  magnitudes  of  the  B-fields  are  1.0  T. 

Fig.  3.  MCD  spectra  with  (a)  normal  magnetic  field  alignment  and  (b)  reversed  magnetic 
field  alignment.  1-at.om  vacancy,  SO  =  17  cm-1,  and  Tmag  =  10  K.  The  magnitudes  of  the 
B-fields  are  1.0  T. 
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Simulated  Absorption  Spectra 


Fig.  la 


Simulated  MCD  Spectra 
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Fig.  lb 
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Normal  Orientation  Magnetic  Field 


Reversed  Magnetic  Field 


Fig.  2a 


Fig.  2b 
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Fig.  3b. 


This  page  is  Distribution  A:  Approved  for  public  release,  distribution  unlimited 

122 


